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Introduction 

An accurate numerical analysis of the flows associated 
with rotor-stator configurations in turhomachincry can he 
very helpful in optimizing the performance of turbomachin- 
erv. Such analyses tend to be computationally expensive 
and extremely complex because of the following reasons: 1) 
the flow is inherently unsteady, 2) the geometries involved 
are complicated, 3) the flow periodically transitions between 
laminar and turbulent flow and 4) there is relative motion 
between the stator and rotor rows. However, a clear under- 
standing of the aerodynamic processes associated with tur- 
boinachinery can aid the design process considerably, and 
hence, the rather large computer costs of simulating the 
three-dimensional unsteady flows associated with turboma- 
chinery are completely justified. 

Several calculations of cascade flow already exist in the 
literature. These studies include two and three-dimensional 
calculations using both the Ruler and Navier-Stokes equa- 
tions. References (1 )-(5) constitute a typical cross-section of 
the work done previously and are by no means a complete 
review of earlier efforts. While analyses of flows through 
isolated rows can be used to study many of the fluid me- 
chanical phenomena in turbomachinery, such analyses do 
not yield any information regarding the unsteadiness aris- 
ing out of the interaction of moving and stationary rows 
of airfoils. These interaction effects become increasingly im- 
portant as the distance between successive rows is decreased. 
The experimental results of (6) show that the temporal pres- 
sure fluctuation near the leading edge of the rotor can be as 
much as 72% of the exit dynamic pressure when the axial 
gap is reduced to 15% of the chord length (for the operating 
conditions and geometry chosen). Thus, the need for treat- 
ing the rotor and stator airfoils as a system in cases where 
interaction effects are predominant is obvious. 

From a computational point of view one major difficulty 
in simulating rotor-stator flows arises because of the relative 
motion of the rotor and stator airfoils. A single grid that 
wraps around both the rotor and stator would have to dis- 
tort considerably to accommodate the motion of the rotor 
and may result in inaccurate calculations. For small values 
of the axial gap between the rotor and stator airfoils such 
an approach may even be altogether impractical. The ob- 
vious solution to this problem is to use several grids which 
move relative to each other. Typically one would use a set 
of stationary grids to envelop the stator airfoils and a set 
of moving grids (stationary with respect to the rotor) to 
envelop the rotor airfoils. Information is then transfered be- 


tween the several grids employed with the help of specialized 
boundary conditions. 

Reference (7) presents rotor-stator interaction results ob- 
tained using the Euler equations. The various natural 
boundary conditions such as inlet, discharge, blade surface 
and periodicity boundary conditions that are required for 
rotor-stator calculations are presented and the unsteady flow 
through a fan stage is calculated. However, there are sev- 
eral areas that have not been addressed in (7). namely, 1) 
a general methodology of information transfer between the 
multiple grids employed, 2) viscous effects and 3) the three- 
dimensionality of the flow. Viscous effects can contribute 
significantly towards the unsteady component of the flow 
because of the passage of the second set of airfoils through 
the wakes of the first set. End wall and tip leakage effects 
and the geometry of the airfoils may contribute significantly 
to the three-dimensionality of the flow. Hence, an accurate 
simulation of the flow within turbomachinery requires the 
time-accurate solution of the unsteady Navier-Stokes equa- 
tions in three-dimensions. 

In the multiple grid approach the calculation is performed 
on several grids that are either patched toget her or are over- 
laid. The boundary conditions used to transfer information 
from one grid to another must satisfy several requirements 
before they can be used effectively. Some of these require- 
ments are listed below. The boundary conditions must be 

1) numerically stable, 

2) spatially and temporally accurate, 

3) easily applicable in generalized coordinates, 

4) conservative so that flow discontinuities can move from 
one grid to another without any distortion. The conservative 
property although desirable is not required in the case of 
flows without discontinuities. 

The boundary conditions required to transfer information 
from patch to patch in the pat died -grid approach are devel- 
oped in detail in ( 8 )-( 10). In (N) a conservative patch bound- 
ary condition is developed for first-order accurate explicit 
schemes. Results demonstrating the conservative property 
of the new boundary condition ami the quality of solutions 
possible with patclied-grids are presented. In (9) and (10), 
this boundary condition is extended to work with implicit 
second-order accurate schemes. The modifications to the 
boundary scheme that are required to transfer in formation 
between two patches that are moving relative to each other 
are also developed in (9) and (10). Preliminary results for a 
rotor-stator configuration are presented in (10). 


The patched grid technique as developed in ( ft )-( 10) is 
used to simulate the flow past a two-dimensional rotor-stator 
configuration of an axial turbine in (11). The airfoil ge- 
ometry and flow conditions used are the same as those in 

(6). The unsteady, thin-layer, Navier-Stokes equations are 
solved in a time accurate manner to obtain the unsteady 
flow field associated with this configuration. The numeri- 
cally obtained results are compared with the experimental 
results of (G). A good comparison between theory and ex- 
periment is obtained in the case of time-averaged pressures 
on the rotor and stator. Pressure amplitudes (corresponding 
to the pressure variation in time) were also found to com- 
pare well with experiment thus indicating the validity of the 
computed unsteady component of the flow. 

In the present study the approximation of 
two-dimensionality is removed and three-dimensional airfoil 
geometries are used. In addition the hub, outer casing and 
rotor tip clearance are all included in the calculation. A sys- 
tem of patched and overlaid grids are used to discretize the 
rather complex geometry of the three-dimensional configu- 
ration. An implicit, upwind third-order accurate method is 
used in all the patches. The equations solved are the un- 
steady, thin-layer Navier-Stokes equations in three dimen- 
sions. The turbulence model used is a modification of the 
Baldwin- Lomax model (12) that is developed in (13). 

The presentation will include a brief description of the 
grid generation procedure, the integration method and the 
various boundary conditions used including the patch and 
overlay boundary conditions. Results in the form of hub to 
tip variations of time-averaged pressures and velocity fields 
to show the various vortical structures (horseshoe vortices 
at the hub and tip for the stator) will also be presented. 
Comparisons with the experimental data of (6) will be made. 
A sample of the results obtained are presented in the next 
section. 

Results 

In this section results obtained for the rotor-stator con- 
figuration shown in Fig. 1 are presented. Approximately 
five cycles (a cycle corresponds to the motion of the rotor 
through an angle equal to 2tt /N where N is the number of 
stator or rotor airfoils) were required to eliminate the ini- 
tial transients and establish a solution that was periodic in 
time. The calculation was performed at a constant time 
step value of about 0.04 (this translates into 2000 time- 
steps/cycle). The inlet Mach number used for this calcu- 
lation was 0.07. The rotor velocity is determined from the 
desired flow coefficient (0.78 in this case) and the inlet axial 
velocity (w^). The Reynolds number used for this calcula- 
tion was 100,000/in. 

Midspan Time-Averaged Pressures 

The available experimental data, for an axial gap of 15% , 
includes time-averaged pressures and pressure amplitudes on 
the rotor and stator airfoils at midspan (G). Figure 2 show's 
the predicted and experimental time-averaged pressure co- 


efficient values ( C p ) as a function of the axial distance along 
the stator. The pressure coefficient is defined as 

y t Pn i ff ~ {pi ) m lc 1 
1 P ” 1 2 

2 Pin let & 

where p avg is the static pressure averaged over one cycle. 
(pt)tniei and p,nUt are the average total pressure and density, 
respectively, at midspan at the inlet and u? is the velocity of 
the rotor at midspan. The comparison between theory and 
experiment is good. A small separation bubble was found on 
the trailing edge circle of the stator in the numerical results. 
This is seen as a spatial fluctuation in pressure towards the 
trailing edge of the st ator. Figure 3 shows the midspan t ime- 
averaged C p distribution for the rotor. As in the case of 
the stator, the comparison between theory and experiment 
is good. A small separation bubble was predicted on the 
trailing edge circle of the rotor. The bubble is seen as a 
sharp dip and rise in the pressure curve. 
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2. l ime- averaged pressure distribution on the stator at 
midspan. 
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3. Time-averaged pressure distribution on the rotor at 
midspan. 
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